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Abstract. We study the dynamics of condensation in a misanthrope process 
with nonlinear jump rates and factorized stationary states. For large enough 
density, it is known that such models have a phase separated state, with a non¬ 
zero fraction of the total mass concentrating in a single lattice site. It has been 
established in [B Waclaw and M R Evans, Phys. Rev. Lett., 108(7):070601, 2012] 
for asymmetric dynamics that such processes exhibit explosive condensation, 
where the time to reach the stationary state vanishes with increasing system 
size. This constitutes a spatially extended version of instantaneous gelation which 
has previously been studied only in mean-field coagulation models. We show 
that this phenomenon also occurs for symmetric dynamics in one dimension if 
the non-linearity is strong enough, and we find a coarsening regime where the 
time to stationarity diverges with the system size for weak non-linearity. In 
higher space dimensions explosive condensation is expected to be generic for all 
parameter values. Our results are based on heuristic mean field arguments which 
are confirmed by simulation data. 
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1. Introduction 

The statistical mechanics of non-equilibrium cluster growth is a topic of recurring 
scientific interest due to a wide range of applications in polymer physics, atmospheric 
science, soft matter and astrophysics (see [1] and the references therein). A large 
collection of clusters is studied, each of which is an aggregate of constituent units 
called monomers. Clusters grow as a result of pairwise interactions. Two basic 
growth mechanisms are common in practice: coagulation and particle exchange. In 
coagulating systems, when a cluster of size i interacts with a cluster of size j, a single 
cluster of size i + j is formed. In particle exchange models, when a cluster of size 
i interacts with a cluster of size j, a monomer is exchanged, resulting in a pair of 
clusters with sizes i — \ and j -I- 1. Related is the well-known Becker-Doring model 
of homogeneous nucleation where monomer exchange is indirect: clusters grow or 
decay by absorbing or shedding one monomer at a time (see [2] and the references 
therein). In both coagulation and particle exchange models, the rates of interaction 
generally depend on the sizes, i and j, of the interacting clusters. This dependence 
is encoded in a (usually symmetric) function, K{i,j), of the cluster sizes called the 
interaction kernel. Given K(i,j), one fundamental question of practical importance 
is to understand the rate at which large clusters are formed. 

A common feature of both processes which is relevant to this question is a 
dynamical transition known as gelation. It occurs when the kernel is a sufficiently 
rapidly increasing function of the particle sizes. Gelation was originally discovered as 
a singularity in the mean-field coagulation equations resulting in the divergence of the 
second moment of the cluster size distribution in a finite time, t*. For a full discussion 
see [3] and the references therein. The mean field equations for exchange driven growth 
can exhibit similar behaviour [4] . The meaning of this singularity becomes clear if one 
considers a finite Markov chain [5] containing N particles and studies its behaviour in 
the limit N ^ oo. This was done explicitly for coagulation with the product kernel, 
K{i,j) = ij, in [6]. It was shown that, for fixed N, a finite fraction of the total mass 
in the system is collected into a single ’’super-particle” within a time tN- As N ^ oo, 
tjv tends to a finite value, t*, known as the gelation time and the mass contained in 
the super-particle diverges. Gelation in an infinite system is therefore interpreted as 
the formation of a cluster of arbitrarily large size within a finite time. 

Soon after the discovery of the gelation transition, it was conjectured [7, 8] that for 
some kernels, the gelation time, t*, is zero. The counter-intuitive idea that arbitrarily 
large clusters could form in arbitrarily short time is called instantaneous gelation. 
This effect was proven rigorously for coagulation models in [9]. A similar regime is 
observed for exchange-driven growth [4] and for the Becker-Doring model [10, II]. 
Since it exhibits a singularity at time 0-I-, instantaneous gelation cannot be described 
using mean field kinetic equations without introducing some kind of regularisation. 
Such regularisation could be done by introducing an explicit cut-off into the kinetic 
theory as was done for coagulation in [12] or by heuristically incorporating a finite 
system size, TV, into the mean field equations as was done for exchange-driven growth 
in [4]. In both cases, there is strong evidence that while the gel time indeed decreases 
as the regularisation is removed, this decrease is logarithmically slow. Consequently 
it is expected to be difficult to distinguish between regular finite time gelation and 
instantaneous gelation in practice. 

Until recently, it was widely believed that instantaneous gelation was a 
pathological property of mean-field kinetic equations and could not occur in spatially 
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extended systems. A counter-example has recently been found [13] in a variation of a 
spatially extended exchange-driven growth model of misanthrope type. Misanthrope 
processes [14, 15] are lattice models in which particles jump between neighbouring 
sites at rates which are a function of the number of particles at both the origin and 
destination sites. If one thinks of each site on the lattice as a cluster, the analogy with 
exchange-driven growth becomes clear. Such processes have been extensively studied 
because they can exhibit stationary measures which can be factorised as a product 
of single site measures and calculated analytically. Furthermore, these stationary 
measures can exhibit an interesting transition as the total particle density is varied. 
If the density exceeds a critical value, a finite fraction of all particles concentrates on 
a single lattice site, whereas the bulk of the system is distributed homogeneously at 
the critical density. This phenomenon is referred to as condensation and the site on 
which a fraction of all particles concentrates is called the condensate. The connection 
between condensation and gelation emerges when one considers the dynamics of how 
the condensate is formed in time. The condensate is the analogue of the super-particle 
discussed above. The question of how the condensate forms dynamically cannot 
be answered solely from knowledge of the stationary product measure itself. The 
model studied in [13] was found to generate a condensate in a time which decreases 
logarithmically with the system size, like in the case of instantaneous gelation. They 
called this effect ’’explosive condensation”. 

The dynamics of condensation in spatially homogeneous particle systems has 
already been studied heuristically in [16] and in subsequent work on zero-range 
processes [17, 18, 19, 20, 21] and related models [22, 23]. In these processes the 
jump rate of particles depends only on the occupation number on the departure site, 
and condensation occurs if this is a sufficiently fast decreasing function. There is also 
a significant literature on the dynamics of condensation in spatially heterogeneous 
models (see [24] and references therein) which is not discussed here. Mathematically 
rigorous results in the homogeneous case have first been obtained for stationary 
distributions (see [25, 19, 26] and references therein), and more recently also on the 
stationary dynamics of the condensate (see e.g. [27, 28]), and the formation of the 
condensate from homogeneous initial conditions [29]. The hydrodynamic behaviour 
has been studied heuristically in [30]. In condensing zero-range type models the 
motion of large clusters is the slowest time scale, exchange of particles occurs through 
fluctuations in size, and relaxation time scales diverge with the system size. 

While misanthrope-type processes share the same factorized stationary states as 
zero-range processes, the dynamics of condensation is entirely different [31, 32]. Large 
clusters move on the same time scale as they exchange particles, or merge when they 
meet. This has first been established for the inclusion process on a rigorous level in 
[33] and with heuristic extensions in [34] . This process has been introduced in [35] and 
is basically a linear version of the explosive condensation model, where condensation 
occurs only due to the scaling of a system parameter [36]. This induces a separation 
of time scales, and relaxation times diverge with the system size. Besides applications 
to energy transport [37], the inclusion process can also be interpreted as a multi¬ 
allele version of the Moran model [38]. The model studied here would be a non-linear 
extension of the latter, where the competitive advantage of a group increases super- 
linearly with the group size. Further applications are discussed in [32] and references 
therein. 

The model of [13] was completely asymmetric in the sense that particles could 
only move to the right. A natural question to ask is whether the results on explosive 
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Figure 1. The two different mechanisms of condensation, with coarsening (left) 
for 7 = 2.25 and explosion (middle and right) for 7 = 4.5, are illustrated using 
the second moment of occupation numbers (31). For explosion, nucleation of 
clusters leads to a gradual increase (zoom on right) until one cluster reaches a 
critical size. In contrast, the coarsening regime (left) is not well separated from 
nucleation and both lead to a gradual increase to reach stationarity. Further 
parameter values are L = 128 and p = 4, d = 0.2 (left) and p = I, d = 0.3 (right). 


condensation depend crucially on the asymmetry of the dynamics which facilitates 
mass transport. In this paper we show that this is not the case, and establish a 
detailed picture of the condensation dynamics for a symmetric system on a one¬ 
dimensional lattice with periodic boundary conditions. If the non-linearity of the 
jump rates is strong enough (characterized by a system parameter 7 > 3) explosive 
condensation occurs analogously to the asymmetric system. Local fluctuations of 
the initial conditions are amplified and clusters form in an initial nucleation regime. 
As soon as one of them reaches a critical size explosion occurs, the cluster moves 
fast enough to cover the whole lattice in a vanishing amount of time and grow into 
the condensate. For weaker non-linearity (with system parameter 7 € (2,3)), the 
critical cluster size cannot be reached and several macroscopic clusters form during 
nucleation. These are then moving across the lattice, exchanging mass in a coarsening 
regime very similar to the dynamics in inclusion processes, and the time to reach a 
single condensate diverges with the system size. These two mechanisms are illustrated 
in Figure 1 using the second moment of the occupation numbers (cf. (31)) as a suitable 
observable to characterize the distance to stationarity. The second mechanism does 
not exist for asymmetric dynamics or symmetric dynamics in higher dimensions. Our 
derivation also applies to asymmetric systems and all dimensions, and provides a slight 
simplification of the arguments presented in [13]. 

The paper is organised as follows. In Section 2 we introduce the model and 
summarize known results on the stationary distributions, emphasizing differences to 
zero-range models regarding the stationary particle current or activity. In Section 3 
we present our main results on the different stages and mechanisms of condensation 
dynamics. We derive the critical cluster size as a key concept, and provide heuristic 
arguments supported by simulation data for explosive and non-explosive condensation 
mechanisms. In Section 4 we conclude with a summary and a short discussion of 
expected results in higher space dimensions. 
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2.1. Definition and notation 

The explosive condensation model {ri{t) : t > 0) is a stochastic particle system defined 
on a lattice of L sites, which we fix to be one-dimensional with periodic boundary 
conditions. Configurations are denoted by ?7 = {rix : x G Al) G Xl where 77^ G N is 
the number of particles at site x G A^. This can be arbitrarily large and the full state 
space is given by . 

A particle at site x jumps to site y ^ x with rate 

c{ri,x,y) := p{x,y)r]f.{d + rjy) with parameters 7 >l,d> 0 . ( 1 ) 

The purely spatial part p{x, y) > 0 is given by rates of an irreducible random walk 
on the lattice A, which determines the dynamics of a single particle together with the 
parameter d > 0. In this paper we focus on symmetric nearest neighbour dynamics 
with 

p{^:y) — 2‘^!/,a;+l + 2^y,x — l ^ (2) 

other choices are shortly discussed in Section 4. The interaction part of the rates is 
such that larger occupation numbers on the departure or target site increase the jump 
rate in a super-linear fashion, controlled by the parameter 7 > 1. These rates are a 
slight variation of the ones introduced in [13], which are 

{{r,x+dr -d'*){d+pyr. (3) 

We present our analysis and numerical results for the rates (1) since they have a 
slightly simpler form. For 7 = 1 both rates coincide, and correspond to the previously 
studied inclusion process [33, 34]. They show the same asymptotic behaviour and all 
results, including factorized steady states as well as condensation dynamics, hold in an 
analogous fashion in both models. For reference, if spatial structure is neglected and 
d is small compared to site occupation numbers, the rate equations for this model are 
equivalent to those studied for exchange-driven growth by Ben-Main and Krapivsky 
[4] with the simple product kernel K(i,j) = (f j)"*'. For reasons, explained below, we 
shall be interested in the regime 7 > 2. At the mean field level, this regime should 
exhibit instantaneous gelation according to the analysis of [4]. 

The dynamics of the continuous-time process (r/(t) : t > 0) is then given by the 
usual master equation 

x,y^A 

for all T 7 S X/,, with the shorthand Ptlv] = ^[Vt = v]- We use the notation with 
77!’^ = rjz — dz,y + 5 z,x for the configuration r/ with one particle moved from site y to 
x, and the usual convention that ptiv] = 0 if 77 ^ < 0 for some z G A. 

2.2. Stationary distributions 

Factorized stationary states for this class of misanthrope processes are well established, 
see [15] as well as [31, 32] and references therein. Since we focus on translation invariant 
systems, we have homogeneous factorized distributions 

= TT with marginals i^Jn] =. (5) 

x&Al ^ ' 
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Here the stationary weights are given by 


w{n) = 

k=l 


{k-ip + d 
k'y 


n as n —)■ c» , 


( 6 ) 


where the asymptotic behaviour is the same as for the rates (3) used in [13]. So the 
single-site partition function is 

OO 

m = E w{k)4>^ < OO for all 4> € [0,1] , (7) 

k=0 

since we assume 7 > 1. The fugacity parameter </> controls the average particle density 
in the grand-canonical setting, which is given by the expectation 

OO 

RW ■= iVx)^ = ki^^[k] = (j)d^ logj^!)) . (8) 

fc =0 

This is a monotone increasing function with i?(0) = 0 and maximal value pc '■= 
i?(l) G (0,oo]. In complete analogy to previous results on zero-range processes 
[16, 17, 18, 19] and the model with rates (3) in [13], the process exhibits condensation 
in the thermodynamic limit if the critical density p^ < 00 , which is the case if and only 
if 7 > 2. In that case, the system phase separates into a homogeneous background 
with density pc, and the excess mass (of order {p — Pc)L particles) concentrates on a 
single lattice site, called the condensate. 

Mathematically, this is formulated in terms of the canonical distributions 

wA) , (9) 


which concentrate on configurations X/, jy = G Xy, : Vx = A a fixed 

particle number X > 0 , and the normalization is simply given by a hnite sum 

over all stationary weights. The process is ergodic on the finite set Xl,n, and ■7rL,N 
is the unique stationary distribution. Condensation can then be understood in terms 
of the equivalence of canonical and grand-canonical ensembles in the thermodynamic 
limit with density p > 0 , 


T^L,N —t 


R{(j)) = p< Pc 
vi, P>Pc 


as X, L —>• OO, N/L —>• p . (10) 


This has been established rigorously in [25, 19, 26, 39] and holds in a weak sense, i.e. 
canonical expectations {f}L,N of bounded, local functions / : Xl,n —> K converge to 
the corresponding expectation (f)^ under the grand-canonical product distribution. 
For p < Pc this also holds for unbounded integrable functions / [31]. 

The interpretation is that for p > pc the system exhibits a phase separated 
(condensed) state with a single condensate as explained above. Depending on the 
initial condition, the condensed phase can have larger initial volume which shrinks 
during time evolution, as is the case in the well-studied coarsening process for zero- 
range dynamics [20, 19, 29] or inclusion processes [33, 34]. The homogeneous (or bulk) 
phase of the system is distributed with a truncated version of the critical distribution 
vi , where occupation numbers are bounded by a size of order L due to the presence of 
the condensed phase. This detail will become important later to compute expectations 
in a phase separated state. 
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2.3. Fundamental diagram 

A typical stationary configuration in the explosive condensation model looks almost 
identical to that of extensively studied condensing zero-range processes with rates 
g{n) = 1 + b/n, since the stationary weights w(n) have the same power-law decay. 
However, as has already been observed in [13] for totally asymmetric rates, the 
dynamics of the process is very different. On the level of stationary distributions 
this is signihed by the expected jump rate 

jL,N ■■= {v2{d + v2))l,n ’ 

and its behaviour in the thermodynamic limit 

jip) ■= jL,N ■ ( 12 ) 

N/L^p 

This corresponds to the activity in a symmetric system, and to the stationary current 
in an asymmetric system where j (p) determines the fundamental diagram of the model. 
For simplicity we refer to this quantity simply as current in the following. With the 
equivalence of ensembles, the limiting current can be computed at least for p < pc 
under the factorized grand-canonical distribution choosing (j) according to ( 10 ), 

■= {v2id + v2))^ = iv2)4>id +' (l^) 

For general zero-range processes it is well known that jgc{4>) = 'Pj therefore j{p) 
is simply given by the inverse of ( 8 ) for p < p^. Since the rates g{n) = 1 -|- 6 /n are 
bounded, (10) also implies that j{p) = 1 for all p > pc. The convergence for finite L 
is dominated by the jump rate g{{p — Pc)L) out of the condensate [40] (as shown in 
Figure 2). So the condensed phase in zero-range models does not contribute to the 
current in this model, and it is known that the condensate remains at a fixed position 
and only moves on a very slow timescale [20, 27, 28]. 

For the model studied here, however, the rates are not bounded and it is clear from 
the power-law decay in ( 6 ) and (13) that the current diverges as p ^ pc. This indicates 
that the condensed phase contributes to, and in fact dominates the stationary current 
for super-critical systems. As opposed to zero-range processes, the condensate moves 
with a high rate of order L'^ as is discussed in more detail later. In the condensed 
state of a system with N particles on L sites where N/L > p^, the canonical stationary 
current is therefore asymptotically given by 

jl.n |(p - PerL'^ipI), + . (14) 

The first contribution is due to the condensate, and the second is a bulk contribution 
determined by expectations w.r.t. the truncated critical distribution i>i, and we 
omitted lower order terms involving the parameter d. Note that {r]2)i = oOi but 
in the phase separated state bulk occupation numbers are bounded by the condensate 
size of order (p — Pc)L, and the truncated moment is therefore also (p2)i < (p — Pc)L- 
Thus 

jL,N ^ 2(p - pc)^+^L'^ -I- const. , (15) 

which diverges as L —>■ oo, and since 7 > 2 the current is dominated by the condensed 
phase. The different behaviour for the fundamental diagrams of zero-range and 
explosive condensation models is illustrated in Figure 2. 

The most striking dynamic effect of this model, however, is the possibility of 
explosive condensation, i.e. the condensed stationary state is reached in a time that 
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Figure 2. Canonical current Jl.n (H) ^ function of density p = N/L for 

different system sizes L (full coloured lines) for the explosive condensation model 
with d = 1 and 7 = 5 (left) and for a zero-range process with rates g(n) = l-\-b/n 
with 6 = 4 (right), resulting from exact recursions in the canonical ensemble (see 
e.g. [40]). It is approximated by the grand-canonical current jgc (13) for p < pc 
(full black line) and by (15) for supercritical densities in the explosive model 
(dashed lines). The finite-size corrections in zero-range processes are also given 
by the current in the condensed phase as explained in the text. Note that both 
plots have very different scales. 


decreases with the system size L. This has been established in [13] for totally 
asymmetric dynamics for all 7 > 2. While there is no a-priori reason to expect 
mean field theory to be relevant to dynamics on a one-dimensional lattice, it is 
interesting to remark that the explosive condensation regime, 7 > 2 , corresponds with 
the instantaneous gelation regime at mean-field level. In the following we investigate 
this phenomenon for symmetric mass transport dynamics. 

3. Main results on condensation dynamics 

In contrast to the totally asymmetric system, which exhibits explosive condensation 
for all 7 > 2, the symmetric model only shows explosion for 7 > 3 and the time to 
stationarity diverges for 7 G (2,3). In the latter regime the approach to stationarity 
is dominated by a cluster coarsening dynamics similar to inculsion processes [33, 34]. 
The main difference between the two scenarios is whether the largest cluster is able 
to reach a critical size rric which has to scale sublinearly with the system size L, as 
is discussed in the next subsection. Once it exceeds that size, it will move across 
the lattice in a time vanishing with system size, and the model exhibits explosion. 
Both cases can be distinguished by a very different scaling of the expected time to 
stationarity {Tss)- This can be defined as 

Tss ■= inf {t > 0 ; > {p - Pc)L} , (16) 

i.e. the time for the process to reach a fully condensed configuration with a condensate 
of size {p — Pc)L- 

3.1. Nucleation and critical cluster size 

The role of the initial conditions is not crucial for our result, as long as they are 
homogeneous. We use two such conditions, for the first one we simply place all 
N particles independently on a uniformly chosen site. This leads to a multinomial 
distribution for ry, and fluctuations of the maximal occupation number are of order 
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Figure 3. Amplification of fluctuations through local dynamics. Left: the 
bimodality of the two-site distribution (17) increases with occupation numbers m 
and the parameter ■y (d = 1). This leads to rough configurations with increasing 
fluctuations as shown on the right for 7 = 7 (d = 1, p = 2, L = 2048) at a very 
early time 3.56 X 10“^, just before one of the clusters starts dominating. 


log L. We also confirmed our results numerically for deterministic initial conditions 
where each site contains exactly p particles for integer densities. This has no 
fluctuations and naturally leads to slower dynamics than the first choice. But also in 
this case the fluctuations in the dynamics lead to clusters of size log L within arbitrarily 
small time intervals, since the dynamics is driven by L independent Poisson processes. 
The resulting equilibration dynamics for both initial conditions shows the same scaling 
and only differs in prefactors. 

The initial fluctuations in occupation numbers are amplified by the local dynamics 
and lead to the nucleation of larger clusters. Consider a two-site system with rates (1) 
and a total of m particles, where 772 = m — rji. With (9) the stationary distribution is 
given by 


7r2.m[?7i, w - rji] = 


w(rji)w{m — 771) 


^2,m 


(17) 


Since the weights w have a power-law tail, the distribution is bimodal with the 
boundaries dominating more and more with increasing m and increasing parameter 
7 as shown in Figure 3. So the mass between any two neighbouring sites tends to 
concentrate in one of them, which leads to rough configurations and the formation of 
clusters as dynamically stable objects. 

These clusters then perform fluctuation activated motion, where a cluster at site 
X jumps to a neighbouring site y = x ±1 with rate 


D{m) = m'^{d + {ri'l))/m ~ m 


7-1 


(18) 


Here we have used that the environment outside the cluster is not yet stationary, so 
{Vy) = 0(1) depending only on the particle density. So the cluster loses a particle at 
rate and the factor l/m is the probability that the remaining particles will 

follow and the whole cluster will move to site y. This results from the probability 
that a symmetric random walk started in site 1 hits site m before site 0 , and this walk 
represents the distribution of mass among both sites corresponding to a path from one 
boundary to the other in Figure 3 (left). The walk is symmetric but continuous-time 
with large jump rates for roughly equal mass distribution, which leads to the bimodal 
structure of the stationary distribution. This also implies that the time for a cluster 
to step is dominated by the time it takes for the first particles to move, while all the 
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others follow within the same time scale. Omitting terms involving the parameter d, 
we have to leading order for the corresponding time scale 

m—1 1 pl/2 


E 


I 




(19) 

J 1/m 


which justifies (18). 

Since 7 > 2 jump rates of clusters increase super-linearly in their size, and through 
their symmetric motion clusters gain mass at rate D(m)/m = Therefore the 

time it takes to grow to a size mi 3> log L which is larger than the initial fluctuations, 
is proportional to 


^ dm = —^ 




/logL 


3-7 


3-7 

7-3 


, 7 e (2,3) 
, 7 > 3 


.( 20 ) 


This vanishes as 7 > 3 and the time to grow a cluster of any size is in fact dominated 
by the maximal initial size logL. For 7 € (2,3) the time diverges depending on the 
final size of the cluster mr, and for macroscopic clusters of size 0{L) it scales like 
L3-7. 

In this context, we define the critical cluster size mr 


/D{mc) = ^~1 


as 

L2/(7-1) ^ 


( 21 ) 


By definition, if m 3> rric the cluster visits all L sites of the lattice instantaneously, 
i.e. in a time o(l) as L —>■ 00 . Together with (20) the system then exhibits explosive 
condensation, where equilibration is dominated by the fast motion of a single largest 
cluster. This is only possible if 

me < T 7 > 3 , (22) 

whereas for 7 € (2,3) we have me ^ L and the system cannot exhibit explosion. 
We discuss both cases in detail below. Note that the nucleation of clusters we just 
discussed is dominated by purely local dynamics, the clusters then start to interact 
globally leading to equilibration of the system. 


S.:S. Coarsening dynamics for 7 G (2, 3) 

If 7 < 3 the critical cluster size me ^ L (21) cannot be reached, and instead several 
clusters of macroscopic size m = 0(L) nucleate from the initial condition within a 
diverging time of order . The system reaches a phase separated state where the 
condensed phase is split in several isolated clusters, which move across the lattice 
and exchange mass in the global part of the dynamics. In this case equilibration is 
therefore dominated by a coarsening process of clusters similar to inclusion processes. 
We analyze this in a mean field approach close to the methods in [34] , where we denote 
by m(t) the average cluster size and by n{t) the number of clusters at time t, such 
that by conservation of mass 

m{t)n{t) = [p-Pe)L . (23) 

Since the background in the phase separated state is stationary, we have to use the 
truncation {r]y)i < m in analogy with (15), to get for the jump rate of a cluster 

D{m) = m'^{d + (rjy))/m ~ mC . 


(24) 
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Figure 4. The expected time to stationarity {Tss) (16) for 7 G (2,3) increases 
with the system size, shown for different values of 7 and d. The predicted scaling 
law (30) with fitted constants is denoted by dashed lines, data are averaged over 
100 realizations. 


Note that this is larger than in the nucleation regime (18) due to enhanced fluctuations 
of the bulk. The average distance between clusters is 

s{t) = L/n{t) = m{t)/{p - Pc) , (25) 

so they meet at a rate proportional to 

D{m{t))/s(t)‘^ = • (26) 


Since the dynamics is symmetric, when two clusters of size mi and m 2 meet they 
merge with probability proportional to l/(mi + m 2 ), in analogy with results for the 
inclusion process [33, 34]. The clusters start exchanging particles as soon as they are 
separated by only one site, and the occupation number of this intermediate site is 
a symmetric random walk in continuous time on the state space {0,... ,mi + m 2 }. 
If the clusters do not merge they still exchange a certain number of particles, which 
happens more frequently and leads to a mass exchange on the same order as merge 
events. This is explained in detail in [34] for inclusion processes which correspond to 
7 = 1 , and applies similarly in the model studied here. 

If two clusters merge, the average size increases by 


Am{t) 


{p- Pc)L _ (p - Pc)L 

n{t) — 1 n(t) 


m{t) 
n(t) — 1 


m{t) . 


(27) 


Putting (26), (27) and the merge probability of order l/m{t) together, the change in 
cluster size is given by 

Ern{t) = const. m{t)'^~‘^ . (28) 

The solution with initial condition m(0) ^ L is given by 


m(t) = C{p 


Pc)(t + io)^^^^ with = -r 

C{p - Pc) 


(29) 


We pull out a factor (p—pc) from the multiplicative constant since simulations support 
the intuition that n{t) is independent of the density, and m{t) scales as p — pc- 

From the scaling law (29) it follows that the time to reach the stationary value 
(p — Pc)L during coarsening scales like . This is of the same order as the time 
scale (20) to grow clusters of macroscopic size, so nucleation and coarsening are not 
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Figure 5. The second moment cr'^{t) increases asymptotically as a power law 
predicted in (32) given by a full line, before it saturates at the stationary value 
(p— pc)^L indicated by dashed lines. Parameter values are 7 = 2.5, d = 0.1, p = 2 
(left), and 7 = 2.25, d = 0.01, p = 3 (right). Data points are averaged over 100 
realizations and errors are comparable to the size of the symbols. Fit constants 
are £0 = 33.6, C = 0.0168 (left) and io = 31.9 and C = 0.545(right), where C 
includes the factor (p — pc)^ in (32). 


well separated in this model (cf. Figure 1), in contrast to previous results for inclusion 
processes where separation results from the scaling of the system parameter d [33, 34]. 
In total, the time to stationarity scales like 

(Tss) ~ ^00 for 7 G (2,3) , (30) 

which is confirmed by simulation data presented in Figure 4. 

To compare the scaling law (29) to numerical data, we use the second moment as 
an observable which is numerically better accessible than the mean cluster size. In a 
phase separated state with n{t) clusters of typical size m{t) this is given by 

'= (z If^^ 

X 

Again, ( 7^)1 = 00 for the critical product distribution with 7 G (2,3), but using the 
truncation by to = 0{L) we see that this term is of order L^~'^ L. So the bulk 
contribution can be neglected for large L, and with (23) we get 

= ip- Pc)m{t) = C{p - pcf{t + . (32) 

This is plotted in Figure 5 and shows good agreement with simulation data in a 
coarsening window, before it saturates at its stationary value {p — Pc)^L. 

3.3. Explosive condensation for 7 > 3 

For 7 > 3 the critical cluster size TOc (21) scales sublinearly with L and 

can be reached by the largest cluster in vanishing time of order (log as discussed 
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Figure 6. The expected time to stationarity (T 55 } (16) for 7 > 3 decreases with 
the system size, shown for p = 2 and different values 7 = 5, d = 0.1 and 1 (left), 
and 7 = 7, d = 0.1 and 1 (right). The predicted scaling law (33) with a fitted 
constant is shown by blue dashed lines, data are averaged over 200 realizations 
and errors are of the size of the symbols. 


in (20). From then on the largest cluster dominates the dynamics, visits all sites of 
the lattice and grows to a condensate containing all excess mass in a time vanishing 
with the system size L. As the speed of the clusters increases with increasing size, the 
time to stationarity is dominated by the nucleation time scale and we expect 

(Tss) ~ (log ^0 for 7 > 3 . (33) 

This is plotted in Figure 6 and shows good agreement with simulation data. 

The time evolution of the average cluster size in this regime is dominated by the 
largest cluster, which gains mass at rate D(m)/m as explained above. Since the bulk 
is still non-stationary we can use (18) for the cluster jump rate, and this leads to 

^m(t) = const. . (34) 

at 

This turns out to coincide with (28) for 7 < 3, even though the dynamics in both 
regimes look very different. For 7 > 3 the right-hand side grows super-linearly in m 
which leads to a finite-time blowup of the solution, 

m{t) =C{p-pc){tbu-t)~^^^'^~^^ with tbu = -^^j^ -r- (35) 

C{p- Pc) 

The blow-up time tbu is determined by the initial condition m(0), and since 77 i( 0 ) ^ 
log L we see that 

tbu ^ (logT)^”^ —^0 as L — 7 . 00 . (36) 

This is consistent with explosive condensation and the scaling of {Tss) in (33). Using 
again the connection cr^(t) = {p — pc)m{t) analogous to (32), we confirm (35) by 
comparison with numerical data in Figure 7. 

The explosive dynamics in this regime is completely analogous to the totally 
asymmetric case studied in [13]. Our derivation of the scaling of {Tss) can be directly 
adapted to that case, and provides an alternative and slightly simpler version of the 
arguments presented there. 
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Figure 7. The second moment cr^{t) follows the prediction from (35) indicated 
by a full gold line for L = 512 with corresponding fitted blow up time given by a 
dash dotted line (left). It saturates at the stationary value (p — pc)^L indicated 
by a dashed red line for L = 32. Plotting against (ti,u — t) in double-log scale 
confirms the functional behaviour following from (35) indicated by straight lines 
(right). Parameter values are 7 = 5, p = 2 and d = 0.1. Data points are averaged 
according to (37) over 50 realizations and errors are of the size of the symbols. 


Note that in contrast to the case 7 < 3, the nucleation and explosive condensation 
regime are well separated, and a transition occurs when the largest cluster reaches 
the critical size as can also be seen in Figure 1. Since the nucleation regime 
is dominated by local dynamics, the time to reach TOc is random and as discussed in 
Section 3.1 is basically determined by a maximum of L Poisson random variables. The 
fluctuations of the latter are on the same scale logL as the expectation, leading to 
persistent randomness even with increasing system size L. Once me is reached and the 
explosion mechanism takes over, relevant observables a'^{t) and the size of the largest 
cluster m(t) follow an almost deterministic trajectory with very small fluctuations (see 
Figure 8 ). Variations in system trajectories therefore result from a random time shift 
corresponding to completion of the nucleation regime, and after that typical individual 
trajectories J2x Vx{t)^ follow (35) where tbu is replaced by a random time (see blue 
curve in Figure 8 ). This implies that the typical behaviour is very different from the 
ensemble average a^{t) (red dashed curve in Figure 8 ). Instead of averaging at fixed 
times t, it is thus more informative to average the time it takes different trajectories 
to reach a level tr^ of the second moment, i.e. 

t{a^) ■= <^inf {s > 0 : > cr^}^ ■ (37) 

X 

Inverting this provides the behaviour of a typical trajectory of the system which is 
used to average the data in Figure 7. 
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Figure 8. Individual trajectories for the largest cluster size shown in black 
compare well with the prediction (35) shown as a thick blue curve, modulo a 
random time shift. The usual ensemble average (dashed red curve) does not 
coincide with the typical behaviour, and representative averages of these dynamics 
are therefore given by (37), which is used in Figure 7. Parameter values are 
L = 128, d = 0.1, and 7 = 3.5, p = 2 (left), and 7 = 5, p = 3 (right). For larger 
7 , fluctuations of individual trajectories get smaller. 


This behaviour is in contrast to the case 7 < 3, where typical trajectories 
fluctuate around the mean with fluctuations decreasing with the system size, due 
to self-averaging effects arising from global coarsening dynamics. This results in a 
decreasing error with L in Figure 4, whereas it is largely independent of L in Figure 

6 . 

3.4- The boundary case 7 = 8 

For the boundary case 7 = 3 in one dimension, our results (28) and (34) predict a 
linear equation 

^rn{t) = const. m{t) , (38) 

with an exponential solution for the growth of the mean cluster size. Analogously, 
this leads to the prediction for the second moment 

(T^(t) = Cie^^* with constants Ci,C 2 > 0 . (39) 

It reaches the stationary value {p—pc)^L in a time of order logL, which is the expected 
scaling of {Tss) in this case. This behaviour is confirmed in Figure 9 for the system 
with the rates (3) originally studied in [13], confirming that our analysis does not 
depend on the detailed functional form of the rates but only on their asymptotic 
behaviour. 

4. Summary and discussion 

We have established a complete picture of condensation dynamics in a symmetric 
version of the model introduced in [13], which includes an explosive regime for 
parameter values 7 > 3, and a novel coarsening regime for 7 S (2,3). The stages 
of the condensation dynamics can be summarized as follows: 
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Figure 9. The second moment cr‘^{t) increases asymptotically exponentially 
with time as predicted in (39) given by a full line, before it saturates at the 
stationary value (p — pc)^L indicated by dashed lines (left). Parameter values 
are 7 = 3, d = 0.1, p = 2, using the rates (3), and fit constants are Ci = 25, 
C 2 = 0.018. This leads to a logarithmic increase for (T 55 } as shown on the right, 

which fits well the prediction log - 1 - ^ logL given by a full line. Data 

o 1 o 2 

points are averaged over 100 realizations and errors are comparable to the size of 
the symbols unless indicated. 


• Initialization: due to fluctuations in the initial conditions or the dynamics at a 
very early stage, the largest occupation numbers in the system scale as log L with 
the system size L. 

• Nucleation: amplification of these initial fluctuations by local dynamics leads to 

large clusters of the critical system size TOc ~ within a time of order 

(logL)^”'*' (20) for 7 > 3. For 7 G (2,3) macroscopic clusters of size 0{L) are 
established in a time of order ( 20 ). 

The remaining time evolution is then governed by two distinct global dynamics: 

• Coarsening for 7 G (2,3): macroscopic clusters move across the lattice and 

exchange particles or merge on a time scale L^~~^ which overlaps with the 
nucleation regime, leading to {Tss) (30). 

• Explosive condensation for 7 > 3: a single cluster dominates the process, covering 
the lattice in vanishing time and growing to a single condensate in an almost 
deterministic fashion. This is clearly separated from and much faster than the 
nucleation regime, which dominates the time scale {Tss) ~ (logT)^”^ (33). 

For the asymmetric transport model, the threshold 7 > 2 for the onset of explosive 
condensation agrees with the threshold for instantaneous gelation in mean field theory 
of exchange-driven growth [4], whereas for symmetric transport, 7 > 3 is required. 
This is consistent with the expectation that spatial models become mean-field like in 
their behaviour only when the spatial transport mechanism causes sufficiently strong 
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mixing to break spatial correlations. It is clear that symmetric transport will mix 
more slowly than asymmetric transport. Our analysis suggests that this difference 
makes possible the existence of the window 2 < 7 < 3 in which spatial correlations 
in the coarsening process in one dimension remain sufficiently strong to prevent the 
explosive behaviour present at mean-field level. 

Our theoretical arguments are slightly different from previous work on asymmetric 
systems [13, 32], but can be directly adapted to the that situation, which leads to 
simpler explanations of the explosive condensation regime. Our results imply that 
asymmetry of the dynamics is in fact not necessary for explosive condensation, and 
that the phenomenon is simply caused by strong enough non-linearity of the jump 
rates. Our arguments can also be directly generalized to higher dimensions. It is 
known that the expected number of steps for a symmetric random walk to visit all 
sites of a d-dimensional lattice with side length L scales like the volume for d > 2, 
with logarthmic corrections for d = 2 and 3. With (21), the critical cluster size then 
scales like 

rOc for d > 2 and 7 > 2 . (40) 

Therefore, we expect explosive condensation for symmetric systems in higher 
dimensions for all values of 7 > 2 , resulting from the scaling of cover times for 
symmetric random walks being similar to asymmetric ones. The non-explosive 
behaviour for 7 G (2, 3) in one dimension is in fact a special case resulting from the 
scaling of cover times and recurrence of random walks. The transient behaviour of 
higher dimensional random walks also implies that clusters gain mass at an increased 
rate D{m) = which leads to a vanishing nucleation time of order (logL)^”'’' 

for all 7 > 2 as a correction to ( 20 ). This is then the expected scaling for {Tss) in 
all higher dimensions. Numerical data of the quality to test predicted scaling laws in 
higher dimensions are hard to obtain and beyond the scope of this paper. The spatial 
part p(x, y) does not have to be nearest neighbour, our results apply directly as long 
as it is irreducible and spatially homogeneous. A further interesting question in this 
context is the role of inhomogeneities, and in particular whether the model could be 
non-explosive on (random) graphs which have bottlenecks and large cover times for 
random walks. 

In general, the simulation of particle systems with non-linearly increasing jump 
rates is computationally very demanding. Rejection based algorithms such as random 
sequential update are highly inefficient, and the only option is Gillespie’s algorithm 
which generates the sample path of general continuous-time Markov chains [41]. 
A single step of a macroscopic cluster requires of order individual transitions 
of particles due to symmetry of the dynamics. This is a factor L more than for 
asymmetric models covered in previous works, which explains the slightly poorer 
quality of our numerical results. Choosing small parameter values for d decreases 
the critical density and leads to a clearer separation between the condensed and the 
fluid phase. This improves the quality of numerical data and also their agreement 
with our theoretical arguments which are based on this separation, which is evident 
from Figures 4 and 6 . It also allows us to replace the full dynamics of a cluster step by 
an effective single transition in the late stages of the dynamics, which is particularly 
efficient in the regime 7 > 3. Details on this can be found in [42], Chapter 5. 
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